Epigenomic Characterization of Lymphoid Neoplasms.

Lymphoid neoplasms represent a heterogeneous group of disease entities and subtypes with markedly different molecular and clinical features. Beyond genetic alterations, lymphoid tumors also show widespread epigenomic changes. These severely affect the levels and distribution of DNA methylation, histone modifications, chromatin accessibility, and three-dimensional genome interactions. DNA methylation stands out as a tracer of cell identity and memory, as B cell neoplasms show epigenetic imprints of their cellular origin and proliferative history, which can be quantified by an epigenetic mitotic clock. Chromatin-associated marks are informative to uncover altered regulatory regions and transcription factor networks contributing to the development of distinct lymphoid tumors. Tumor-intrinsic epigenetic and genetic aberrations cooperate and interact with microenvironmental cells to shape the transcriptome at different phases of lymphoma evolution, and intraclonal heterogeneity can now be characterized by single-cell profiling. Finally, epigenetics offers multiple clinical applications, including powerful diagnostic and prognostic biomarkers as well as therapeutic targets.

BCL3 rearrangements in B-cell lymphoid neoplasms occur in two breakpoint clusters associated with different diseases.

The t(14;19)(q32;q13) often juxtaposes BCL3 with immunoglobulin heavy chain (IGH) resulting in overexpression of the gene. In contrast to other oncogenic translocations, BCL3 rearrangement (BCL3-R) has been associated with a broad spectrum of lymphoid neoplasms. Here we report an integrative whole-genome sequence, transcriptomic, and DNA methylation analysis of 13 lymphoid neoplasms with BCL3-R. The resolution of the breakpoints at single base-pair revealed that they occur in two clusters at 5' (n=9) and 3' (n=4) regions of BCL3 associated with two different biological and clinical entities. Both breakpoints were mediated by aberrant class switch recombination of the IGH locus. However, the 5' breakpoints (upstream) juxtaposed BCL3 next to an IGH enhancer leading to overexpression of the gene whereas the 3' breakpoints (downstream) positioned BCL3 outside the influence of the IGH and were not associated with its expression. Upstream BCL3-R tumors had unmutated IGHV, trisomy 12, and mutated genes frequently seen in chronic lymphocytic leukemia (CLL) but had an atypical CLL morphology, immunophenotype, DNA methylome, and expression profile that differ from conventional CLL. In contrast, downstream BCL3-R neoplasms were atypical splenic or nodal marginal zone lymphomas (MZL) with mutated IGHV, complex karyotypes and mutated genes typical of MZL. Two of the latter four tumors transformed to a large B-cell lymphoma. We designed a novel fluorescence in situ hybridization assay that recognizes the two different breakpoints and validated these findings in 17 independent tumors. Overall, upstream or downstream breakpoints of BCL3-R are mainly associated with two subtypes of lymphoid neoplasms with different (epi)genomic, expression, and clinicopathological features resembling atypical CLL and MZL, respectively.

MALAT1 expression is associated with aggressive behavior in indolent B-cell neoplasms.

MALAT1 long non-coding RNA has oncogenic roles but has been poorly studied in indolent B-cell neoplasms. Here, MALAT1 expression was analyzed using RNA-seq, microarrays or qRT-PCR in primary samples from clinico-biological subtypes of chronic lymphocytic leukemia (CLL, n = 266), paired Richter transformation (RT, n = 6) and follicular lymphoma (FL, n = 61). In peripheral blood (PB) CLL samples, high MALAT1 expression was associated with a significantly shorter time to treatment independently from other known prognostic factors. Coding genes expressed in association with MALAT1 in CLL were predominantly related to oncogenic pathways stimulated in the lymph node (LN) microenvironment. In RT paired samples, MALAT1 levels were lower, concordant with their acquired increased independency of external signals. Moreover, MALAT1 levels in paired PB/LN CLLs were similar, suggesting that the prognostic value of MALAT1 expression in PB is mirroring expression differences already present in LN. Similarly, high MALAT1 expression in FL predicted for a shorter progression-free survival, in association with expression pathways promoting FL pathogenesis. In summary, MALAT1 expression is related to pathophysiology and more aggressive clinical behavior of indolent B-cell neoplasms. Particularly in CLL, its levels could be a surrogate marker of the microenvironment stimulation and may contribute to refine the clinical management of these patients.

Whole-genome sequencing of chronic lymphocytic leukemia identifies subgroups with distinct biological and clinical features.

The value of genome-wide over targeted driver analyses for predicting clinical outcomes of cancer patients is debated. Here, we report the whole-genome sequencing of 485 chronic lymphocytic leukemia patients enrolled in clinical trials as part of the United Kingdom's 100,000 Genomes Project. We identify an extended catalog of recurrent coding and noncoding genetic mutations that represents a source for future studies and provide the most complete high-resolution map of structural variants, copy number changes and global genome features including telomere length, mutational signatures and genomic complexity. We demonstrate the relationship of these features with clinical outcome and show that integration of 186 distinct recurrent genomic alterations defines five genomic subgroups that associate with response to therapy, refining conventional outcome prediction. While requiring independent validation, our findings highlight the potential of whole-genome sequencing to inform future risk stratification in chronic lymphocytic leukemia.

Detection of early seeding of Richter transformation in chronic lymphocytic leukemia.

Richter transformation (RT) is a paradigmatic evolution of chronic lymphocytic leukemia (CLL) into a very aggressive large B cell lymphoma conferring a dismal prognosis. The mechanisms driving RT remain largely unknown. We characterized the whole genome, epigenome and transcriptome, combined with single-cell DNA/RNA-sequencing analyses and functional experiments, of 19 cases of CLL developing RT. Studying 54 longitudinal samples covering up to 19 years of disease course, we uncovered minute subclones carrying genomic, immunogenetic and transcriptomic features of RT cells already at CLL diagnosis, which were dormant for up to 19 years before transformation. We also identified new driver alterations, discovered a new mutational signature (SBS-RT), recognized an oxidative phosphorylation (OXPHOS)high-B cell receptor (BCR)low-signaling transcriptional axis in RT and showed that OXPHOS inhibition reduces the proliferation of RT cells. These findings demonstrate the early seeding of subclones driving advanced stages of cancer evolution and uncover potential therapeutic targets for RT.

Molecular map of chronic lymphocytic leukemia and its impact on outcome.

Recent advances in cancer characterization have consistently revealed marked heterogeneity, impeding the completion of integrated molecular and clinical maps for each malignancy. Here, we focus on chronic lymphocytic leukemia (CLL), a B cell neoplasm with variable natural history that is conventionally categorized into two subtypes distinguished by extent of somatic mutations in the heavy-chain variable region of immunoglobulin genes (IGHV). To build the 'CLL map,' we integrated genomic, transcriptomic and epigenomic data from 1,148 patients. We identified 202 candidate genetic drivers of CLL (109 new) and refined the characterization of IGHV subtypes, which revealed distinct genomic landscapes and leukemogenic trajectories. Discovery of new gene expression subtypes further subcategorized this neoplasm and proved to be independent prognostic factors. Clinical outcomes were associated with a combination of genetic, epigenetic and gene expression features, further advancing our prognostic paradigm. Overall, this work reveals fresh insights into CLL oncogenesis and prognostication.

IGLV3-21R110 identifies an aggressive biological subtype of chronic lymphocytic leukemia with intermediate epigenetics.

B-cell receptor (BCR) signaling is crucial for chronic lymphocytic leukemia (CLL) biology. IGLV3-21-expressing B cells may acquire a single point mutation (R110) that triggers autonomous BCR signaling, conferring aggressive behavior. Epigenetic studies have defined 3 CLL subtypes based on methylation signatures reminiscent of naïve-like (n-CLL), intermediate (i-CLL), and memory-like (m-CLL) B cells with different biological features. i-CLL carries a borderline IGHV mutational load and significantly higher use of IGHV3-21/IGLV3-21. To determine the clinical and biological features of IGLV3-21R110 CLL and its relationship to these epigenetic subtypes, we characterized the immunoglobulin gene of 584 CLL cases using whole-genome/exome and RNA sequencing. IGLV3-21R110 was detected in 6.5% of cases: 30 (38%) of 79 i-CLLs, 5 (1.7%) of 291 m-CLLs, and 1 (0.5%) of 189 n-CLLs. All stereotype subset 2 cases carried IGLV3-21R110, whereas 62% of IGLV3-21R110 i-CLL cases had nonstereotyped BCR immunoglobulins. IGLV3-21R110 i-CLL had a significantly higher number of SF3B1 and ATM mutations and total number of driver alterations. However, the R110 mutation was the sole alteration in 1 i-CLL and was accompanied only by del(13q) in 3. Although IGHV mutational status varied, IGLV3-21R110 i-CLL transcriptomically resembled n-CLL/unmutated IGHV CLL with a specific signature including WNT5A/B overexpression. In contrast, i-CLL lacking IGLV3-21R110 mirrored m-CLL/mutated IGHV. Patients with IGLV3-21R110 i-CLL had a short time to first treatment and overall survival similar to those of n-CLL/unmutated IGHV patients, whereas patients with non-IGLV3-21R110 i-CLL had a good prognosis similar to that of patients with m-CLL/mutated IGHV. IGLV3-21R110 defines a CLL subgroup with specific biological features and an unfavorable prognosis independent of IGHV mutational status and epigenetic subtype.

Chromatin activation as a unifying principle underlying pathogenic mechanisms in multiple myeloma.

Multiple myeloma (MM) is a plasma cell neoplasm associated with a broad variety of genetic lesions. In spite of this genetic heterogeneity, MMs share a characteristic malignant phenotype whose underlying molecular basis remains poorly characterized. In the present study, we examined plasma cells from MM using a multi-epigenomics approach and demonstrated that, when compared to normal B cells, malignant plasma cells showed an extensive activation of regulatory elements, in part affecting coregulated adjacent genes. Among target genes up-regulated by this process, we found members of the NOTCH, NF-kB, MTOR signaling, and TP53 signaling pathways. Other activated genes included sets involved in osteoblast differentiation and response to oxidative stress, all of which have been shown to be associated with the MM phenotype and clinical behavior. We functionally characterized MM-specific active distant enhancers controlling the expression of thioredoxin (TXN), a major regulator of cellular redox status and, in addition, identified PRDM5 as a novel essential gene for MM. Collectively, our data indicate that aberrant chromatin activation is a unifying feature underlying the malignant plasma cell phenotype.

Genomic and epigenomic insights into the origin, pathogenesis, and clinical behavior of mantle cell lymphoma subtypes.

Mantle cell lymphoma (MCL) is a mature B-cell neoplasm initially driven by CCND1 rearrangement with 2 molecular subtypes, conventional MCL (cMCL) and leukemic non-nodal MCL (nnMCL), that differ in their clinicobiological behavior. To identify the genetic and epigenetic alterations determining this diversity, we used whole-genome (n = 61) and exome (n = 21) sequencing (74% cMCL, 26% nnMCL) combined with transcriptome and DNA methylation profiles in the context of 5 MCL reference epigenomes. We identified that open and active chromatin at the major translocation cluster locus might facilitate the t(11;14)(q13;32), which modifies the 3-dimensional structure of the involved regions. This translocation is mainly acquired in precursor B cells mediated by recombination-activating genes in both MCL subtypes, whereas in 8% of cases the translocation occurs in mature B cells mediated by activation-induced cytidine deaminase. We identified novel recurrent MCL drivers, including CDKN1B, SAMHD1, BCOR, SYNE1, HNRNPH1, SMARCB1, and DAZAP1. Complex structural alterations emerge as a relevant early oncogenic mechanism in MCL, targeting key driver genes. Breakage-fusion-bridge cycles and translocations activated oncogenes (BMI1, MIR17HG, TERT, MYC, and MYCN), generating gene amplifications and remodeling regulatory regions. cMCL carried significant higher numbers of structural variants, copy number alterations, and driver changes than nnMCL, with exclusive alterations of ATM in cMCL, whereas TP53 and TERT alterations were slightly enriched in nnMCL. Several drivers had prognostic impact, but only TP53 and MYC aberrations added value independently of genomic complexity. An increasing genomic complexity, together with the presence of breakage-fusion-bridge cycles and high DNA methylation changes related to the proliferative cell history, defines patients with different clinical evolution.

The proliferative history shapes the DNA methylome of B-cell tumors and predicts clinical outcome.

We report a systematic analysis of the DNA methylation variability in 1,595 samples of normal cell subpopulations and 14 tumor subtypes spanning the entire human B-cell lineage. Differential methylation among tumor entities relates to differences in cellular origin and to de novo epigenetic alterations, which allowed us to build an accurate machine learning-based diagnostic algorithm. We identify extensive patient-specific methylation variability in silenced chromatin associated with the proliferative history of normal and neoplastic B cells. Mitotic activity generally leaves both hyper- and hypomethylation imprints, but some B-cell neoplasms preferentially gain or lose DNA methylation. Subsequently, we construct a DNA methylation-based mitotic clock called epiCMIT, whose lapse magnitude represents a strong independent prognostic variable in B-cell tumors and is associated with particular driver genetic alterations. Our findings reveal DNA methylation as a holistic tracer of B-cell tumor developmental history, with implications in the differential diagnosis and prediction of clinical outcome.

DNA methylation profiles in chronic lymphocytic leukemia patients treated with chemoimmunotherapy.

BACKGROUND: In order to gain insight into the contribution of DNA methylation to disease progression of chronic lymphocytic leukemia (CLL), using 450K Illumina arrays, we determined the DNA methylation profiles in paired pre-treatment/relapse samples from 34 CLL patients treated with chemoimmunotherapy, mostly (n = 31) with the fludarabine-cyclophosphamide-rituximab (FCR) regimen. RESULTS: The extent of identified changes in CLL cells versus memory B cells from healthy donors was termed "epigenetic burden" (EB) whereas the number of changes between the pre-treatment versus the relapse sample was termed "relapse changes" (RC). Significant (p < 0.05) associations were identified between (i) high EB and short time-to-first-treatment (TTFT); and, (ii) few RCs and short time-to-relapse. Both the EB and the RC clustered in specific genomic regions and chromatin states, including regulatory regions containing binding sites of transcription factors implicated in B cell and CLL biology. CONCLUSIONS: Overall, we show that DNA methylation in CLL follows different dynamics in response to chemoimmunotherapy. These epigenetic alterations were linked with specific clinical and biological features.

Integrated epigenomic and transcriptomic analysis reveals TP63 as a novel player in clinically aggressive chronic lymphocytic leukemia.

Chronic lymphocytic leukemia (CLL) stereotyped subsets #6 and #8 include cases expressing unmutated B cell receptor immunoglobulin (BcR IG) (U-CLL). Subset #6 (IGHV1-69/IGKV3-20) is less aggressive compared to subset #8 (IGHV4-39/IGKV1(D)-39) which has the highest risk for Richter's transformation among all CLL. The underlying reasons for this divergent clinical behavior are not fully elucidated. To gain insight into this issue, here we focused on epigenomic signatures and their links with gene expression, particularly investigating genome-wide DNA methylation profiles in subsets #6 and #8 as well as other U-CLL cases not expressing stereotyped BcR IG. We found that subset #8 showed a distinctive DNA methylation profile compared to all other U-CLL cases, including subset #6. Integrated analysis of DNA methylation and gene expression revealed significant correlation for several genes, particularly highlighting a relevant role for the TP63 gene which was hypomethylated and overexpressed in subset #8. This observation was validated by quantitative PCR, which also revealed TP63 mRNA overexpression in additional nonsubset U-CLL cases. BcR stimulation had distinct effects on p63 protein expression, particularly leading to induction in subset #8, accompanied by increased CLL cell survival. This pro-survival effect was also supported by siRNA-mediated downregulation of p63 expression resulting in increased apoptosis. In conclusion, we report that DNA methylation profiles may vary even among CLL patients with similar somatic hypermutation status, supporting a compartmentalized approach to dissecting CLL biology. Furthermore, we highlight p63 as a novel prosurvival factor in CLL, thus identifying another piece of the complex puzzle of clinical aggressiveness.

The reference epigenome and regulatory chromatin landscape of chronic lymphocytic leukemia.

Chronic lymphocytic leukemia (CLL) is a frequent hematological neoplasm in which underlying epigenetic alterations are only partially understood. Here, we analyze the reference epigenome of seven primary CLLs and the regulatory chromatin landscape of 107 primary cases in the context of normal B cell differentiation. We identify that the CLL chromatin landscape is largely influenced by distinct dynamics during normal B cell maturation. Beyond this, we define extensive catalogues of regulatory elements de novo reprogrammed in CLL as a whole and in its major clinico-biological subtypes classified by IGHV somatic hypermutation levels. We uncover that IGHV-unmutated CLLs harbor more active and open chromatin than IGHV-mutated cases. Furthermore, we show that de novo active regions in CLL are enriched for NFAT, FOX and TCF/LEF transcription factor family binding sites. Although most genetic alterations are not associated with consistent epigenetic profiles, CLLs with MYD88 mutations and trisomy 12 show distinct chromatin configurations. Furthermore, we observe that non-coding mutations in IGHV-mutated CLLs are enriched in H3K27ac-associated regulatory elements outside accessible chromatin. Overall, this study provides an integrative portrait of the CLL epigenome, identifies extensive networks of altered regulatory elements and sheds light on the relationship between the genetic and epigenetic architecture of the disease.

Decoding the DNA Methylome of Mantle Cell Lymphoma in the Light of the Entire B Cell Lineage.

We analyzed the in silico purified DNA methylation signatures of 82 mantle cell lymphomas (MCL) in comparison with cell subpopulations spanning the entire B cell lineage. We identified two MCL subgroups, respectively carrying epigenetic imprints of germinal-center-inexperienced and germinal-center-experienced B cells, and we found that DNA methylation profiles during lymphomagenesis are largely influenced by the methylation dynamics in normal B cells. An integrative epigenomic approach revealed 10,504 differentially methylated regions in regulatory elements marked by H3K27ac in MCL primary cases, including a distant enhancer showing de novo looping to the MCL oncogene SOX11. Finally, we observed that the magnitude of DNA methylation changes per case is highly variable and serves as an independent prognostic factor for MCL outcome.

Whole-genome fingerprint of the DNA methylome during human B cell differentiation.

We analyzed the DNA methylome of ten subpopulations spanning the entire B cell differentiation program by whole-genome bisulfite sequencing and high-density microarrays. We observed that non-CpG methylation disappeared upon B cell commitment, whereas CpG methylation changed extensively during B cell maturation, showing an accumulative pattern and affecting around 30% of all measured CpG sites. Early differentiation stages mainly displayed enhancer demethylation, which was associated with upregulation of key B cell transcription factors and affected multiple genes involved in B cell biology. Late differentiation stages, in contrast, showed extensive demethylation of heterochromatin and methylation gain at Polycomb-repressed areas, and genes with apparent functional impact in B cells were not affected. This signature, which has previously been linked to aging and cancer, was particularly widespread in mature cells with an extended lifespan. Comparing B cell neoplasms with their normal counterparts, we determined that they frequently acquire methylation changes in regions already undergoing dynamic methylation during normal B cell differentiation.